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We investigate topological properties of density matrices motivated by the question to what extent phenomena 
such as topological insulators and superconductors can be generalized to mixed states in the framework of open 
quantum systems. The notion of geometric phases has been extended from pure to mixed states by Uhlmann in 
[Rep. Math. Phys. 24, 229 (1986)], where an emergent gauge theory over the density matrices based on their 
pure-state representation in a larger Hilbert space has been reported. However, since the uniquely defined square 
root yof a density matrix p provides a global gauge, this construction is always topologically trivial. Here, 
we study a more restrictive gauge structure which can be topologically non-trivial and is capable of resolving 
homotopically distinct mappings of density matrices subject to various spectral constraints. Remarkably, in 
this framework, topological invariants can be directly defined and calculated for mixed states. In the limit of 
pure states, the well known system of topological invariants for gapped band structures at zero temperature is 
reproduced. We compare our construction with recent approaches to Chem insulators at finite temperature. 


I. INTRODUCTION 

With the identification of geometric phases in quantum 
physics [1-3], an emergent classical gauge structure in 
elementary quantum mechanics has been revealed. Several 
fundamental discoveries of contemporary quantum physics 
such as the integer quantum Hall effect [4-7] and, more 
recently, the periodic table of topological insulators and 
superconductors [8-12] can be theoretically described as 
topological properties of such emergent gauge theories in the 
framework of Bloch bands: The local U(n) gauge degree 
of freedom there consists in the choice of an arbitrary or¬ 
thonormal basis of Bloch states \u%) that span the projection 
P(k) = Y^a=i\ u k)( u k I onto t ^ le n occupied bands of an in¬ 
sulating band structure at lattice momentum k. The covariant 
derivative for such gauge theories has been constructed by 
Kato as early as 1950 in his proof of the adiabatic theorem 
of quantum mechanics [13], Topological invariants of this 
gauge structure distinguish homotopically distinct mappings 
k i—^ P(fc), i.e., topologically inequivalent band structures. 
The consideration of generic symmetries [14] refines this 
homotopy classification and leads to the periodic table of 
topological invariants [8-12] (see Ref. [15] for a recent 
rigorous discussion). 

Real physical systems, however, are not perfectly isolated 
from their environment and are to be described by a mixed 
state density matrix p rather than a pure state wave function. 
It is hence natural to ask whether the concept of geometric 
phases in quantum physics can be generalized to the realm 
of density matrices and their quantum mechanical evolution. 
This question has been answered in the affirmative in a series 
of papers by Uhlmann [16-18], who viewed the redundancy 
in the purification of density matrices, i.e., in the representa¬ 
tion of p as a pure state in a larger Hilbert space, as a gauge 
degree of freedom: p = ww* = (wU)(wliy, where the 
Hilbert Schmidt operator w denotes a purification of p and U 
is a unitary operator. In this context, the Uhlmann-connection 


defining parallel transport of Hilbert Schmidt operators along 
a path t i—^ p(t) of density matrices is given by the constraint 
[16] 

w'w — w'w = 0, (1) 

where the dot denotes the derivative with respect to t. 

The purpose of our present work is to unravel the topolog¬ 
ical properties of various gauge structures over the space of 
quantum mechanical density matrices. Since the purification 
w = ^J~p, uniquely defined for every density matrix p, pro¬ 
vides a global gauge, the general purification scheme in terms 
of Hilbert Schmidt operators as considered by Uhlmann is 
always topologically trivial. However, we would like to point 
out that this scheme can also be applied to density matrices of 
pure states which then would also yield a topologically trivial 
gauge theory - despite the fact that pure states give rise to a 
rich spectrum of topological phenomena as mentioned above. 
Thus, the fact that some topologically trivial gauge structure 
can be constructed over the space of density matrices by no 
means implies that there are no topological features to be 
discovered. In fact, more generally speaking, every topo¬ 
logical gauge theory can be embedded into a topologically 
trivial theory with a larger gauge degree of freedom over 
the same base space [19], This raises the natural question 
whether a more restricted notion than general purification 
of density matrices in terms of Hilbert Schmidt operators 
could be employed to reveal topological aspects of families of 
density matrices p(k) parameterized by a lattice momentum k. 

Key results - We investigate an ensemble of pure states 
(EOPS) scheme in terms of non-orthonormalized pure states 
\ip a ) satisfying 

p = £ irxri (2) 

a 

as a gauge structure over the space of density matrices p. This 
construction draws intuition from the entanglement theory of 
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mixed states where, e.g., the entanglement of formation [ 20 ] 
is defined in terms of similar ensembles of pure states. Start¬ 
ing from the spectral representation p = '^2 a p a \'it> a )(ip a \> 
a natural EOPS is \ip a ) = y/p^\ij) a ). Under arbitrary uni¬ 
tary rotations of this frame, the pure states are no longer mu¬ 
tually orthogonal but still project down to the same density 
matrix p via Eq. (2). The resulting gauge theory of non- 
orthogonal frames over the density matrices is a generaliza¬ 
tion of the pure state case. For invertible density matrices 
without any further spectral assumptions, it is equivalent to 
Uhlmann’s construction and thus topologically trivial. Inter¬ 
estingly, if assumptions on the spectral degeneracy of the den¬ 
sity matrices are made, the present scheme can accommodate 
topologically non-trivial mixed states. Topologically inequiv¬ 
alent mappings k K>■ p(k) in this framework cannot be contin¬ 
uously deformed into each other without either violating the 
spectral assumptions or breaking possible protecting physical 
symmetries. Our construction formalizes the notion of a pu¬ 
rity gap proposed in Refs. [21, 22] - a purity gap closing 
is a level-crossing in p and thus a violation of the underlying 
spectral assumptions that makes a change in the topology pos¬ 
sible. Under non-equilibrium conditions, a physical system 
may be non-ergodic and may thus be described by a singular 
density matrix of rank n < N, where N is the dimension of 
the Hilbert space. Also in this case, the EOPS gauge struc¬ 
ture can become topologically non-trivial. The topological 
invariants associated with the proposed gauge structure are 
constructed in terms of Uhlmann’s connection (see Eq. (1)) 
without reference to the pure state case which is, however, 
correctly reproduced as a limiting case here. 

Finally, we discuss notable differences to recent publica¬ 
tions [23-25] where topological invariants for two-banded 
one-dimensional systems and Chern insulators [26] at finite 
temperature, respectively, have been proposed that are not 
characteristic classes of a gauge theory. We illustrate with 
an explicit example that a two-dimensional system at finite 
temperature is in general not uniquely characterized by a 
single invariant of that kind. 

Outline- This article is structured as follows. In Section 
II, we review how a gauge structure emerges for both pure 
and mixed states in elementary quantum mechanics. In Sec¬ 
tion III, topological invariants for density matrices with vari¬ 
ous spectral assumptions are defined. We compare the present 
approach to the construction of Uhlmann phase winding num¬ 
bers that have recently been proposed [24, 25] to classify ther¬ 
mal Chern insulators in Section IV. Concluding remarks are 
presented in Section V. 


H. EMERGENT GAUGE STRUCTURES FROM PURE TO 
MIXED STATES 

A. Adiabatic time evolution of pure states 

A natural and conceptually simple scenario for the occur¬ 
rence of geometric phases is the time evolution of a quan¬ 
tum mechanical system with a Hamiltonian H(R(t)) that de¬ 


pends adiabatically on time via some control parameters R(t) 
[27], As a side remark, we note that the concept of geometric 
phases has been generalized to non-adiabatic [28] and non- 
cyclic [29, 30] evolution. The non-degenerate ground state of 
H(R) is denoted by | R). Since the system cannot leave its 
instantaneous ground state when initially prepared in |f?(fo)), 
the geometric constraint 

P{t)Mt)) = |*(*)> ( 3 ) 

with P(t) = |f?(t))(i?(f)| the projection onto the ground 
state, is imposed on the solution |’&(t)) of the time dependent 
Schrodinger equation in the adiabatic limit. This constraint, 
reflecting the absence of any dynamical level transitions, im¬ 
plies the form |4'(f)) = e*^W - ^ D W^|i?(f)), where \R(t)) is 
a family of ground states with a smooth relative phase, i.e., a 
gauge of ground states. </>£>(() = drE q(t) with the ground 
state energy Eq is called the dynamical phase. The additional 
phase factor of the state vector \i/t(t)) = e*^ D ( t )|'I'(i)) 
relative to the gauge \R(t)) reveals a deep geometric “prin¬ 
ciple of least effort” that nature employs in the adiabatic 
limit of quantum mechanics: \ip(t)) is the shortest path in 
Hilbert space that satisfies Eq. (3), shortest as measured 

by L = dry (i/>(t)|t/>(t)) with “velocity” \ip (r)), i.e., 
by the metric induced by the inner product in Hilbert space. 
To systematically construct this path, it is helpful to decom¬ 
pose its tangent vector into two orthogonal components as 
\ip{t)) = P(t)\ip(t)) + (1 — P(t))\ip(t)). The first compo¬ 
nent is called the vertical part \ip)v and generates an evolu¬ 
tion that stays within the projection P(t), i.e., a mere change 
of phase of |^>(f)). The second part in contrast generates an 
evolution perpendicular to P(t) in the sense of the inner prod¬ 
uct and is consequently called the horizontal part \iP)h- From 

IIV’II 2 = HVvIl 2 + HV’irll 2 , it is easy to see that the tangent 
vector to the shortest path must be purely horizontal, i.e., 

P(fM*)}=0, (4) 

or, equivalently for non-degenerate ground states, 

(ij}(t)\ijj(t)) = 0. This parallel-transport prescription 

immediately determines 

<t>(t) = i dT ( R ( T )\^ T \R('r)). (5) 

The quantity Ad^(R(t)) = (i?(t)|4:|i?(f)) transforms un- 
der a gauge transformation \R) — > e**^! R) as Ad_ —> 
Ajl + i-ffix(R(f))’ i- e -’ as a y(l) g au g e field. The geomet¬ 
ric phase </> 7 = i A (mod 27r) associated with a loop 7 in 
parameter space is gauge invariant. This construction can be 
immediately generalized to n-fold degenerate ground states 
where the projection P(t) in Eq. (4) becomes 

n 

pw=mp(t)=Y,\ Rj m Rj (t)\ ( 6 ) 

3 =1 

with the ground state manifold basis or frame / = 

(ji? 1 ),..., | R n )). We note the independence of the projection 
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under U(n) basis transformations / —► fU. Thus, instead of 
a mere phase factor, the local gauge degree of freedom is then 
a U ( n ) basis transformation on the ground state manifold. In 
this case, the geometric phase or, in mathematical terms, the 
(7(1) holonomy e u ^ is replaced by the U{n) holonomy 

U'i = Te~ f i A , (7) 

where A-* 1 = (B?\d\R l ) is the non-Abelian gauge field or con¬ 
nection and T is the time ordering operator along the cyclic 
path 7. 

We finally note that the connection can be expressed in 
a manifestly gauge invariant or basis independent form as 
Ak = — [(c7P),P] [13]. In this case, the geometric phase 
in the ground state manifold for a loop 7 starting at time to 
is described by the propagator of the Schrodinger equation in 
the adiabatic limit, or the holonomy operator 

Ul = P{to)Te~^ AK P{to). ( 8 ) 

U 1 as occurring in Eq. (7) is then just the gauge (basis) 
dependent representation matrix of U]^. 


B. From geometric phases to topological band structures 

Let us consider an insulating band structure in d spatial di¬ 
mensions with n occupied Bloch states |u£) below the Fermi 
energy that span the projection P(k) = Y^,=i\ u k) ( u k\ aflat - 
tice momentum k. If we identify the energy gap relevant for 
the adiabatic approximation with the band gap of the insulator 
and the parameter manifold with the Brillouin zone (BZ) in 
which the lattice momentum is defined, a gauge structure can 
be defined in complete analogy to Section II A. The projection 
P(k) defines an n- plane in the A'-dimensional Hilbert space 
defined at every momentum k in the BZ, where N is the total 
number of bands (occupied and empty) of the model system. 
P(k) is thus by definition a point on the Grassmann mani¬ 
fold G n (C N ) = U{N)/{U{n) x U{N — n)). Topological 
invariants of this gauge structure then distinguish topologi¬ 
cally inequivalent mappings k i-> P{k), that is homotopically 
inequivalent mappings from the BZ torus T d to G n (C N ). For 
the case of particle number conserving insulators without ad¬ 
ditional symmetries outlined here, an integer invariant named 
the m-th Chern number [31] is defined for even spatial dimen¬ 
sions d = 2m while all insulators are equivalent in odd spa¬ 
tial dimensions. Physical symmetries imply constraints on the 
form of the mapping k 1 —>• P(k). For example the anti-unitary 
time reversal symmetry T yields P(k) = TP{— fc)T -1 . Map¬ 
pings (band structures) that would be topologically equiva¬ 
lent when breaking the relevant symmetries may be topolog¬ 
ically inequivalent when maintaining the symmetries. This 
phenomenon is commonly referred to as symmetry protection 
of a topological state. Taking also into account superconduct¬ 
ing band structures as well as all symmetries considered by 
Altland and Zirnbauer [14] significantly refines the system of 
topological invariants and leads to a pattern which has been 


coined the periodic table of topological insulators and super¬ 
conductors [8-12], In this work, we would like to address the 
question to what extent such invariants can be generalized to 
lattice translation-invariant quantum many body systems in a 
mixed state. While the relevant physical symmetries can read¬ 
ily be generalized to the realm of density matrices, the defini¬ 
tion of topological invariants for mixed states from a gauge 
structure over the space of density matrices has not been dis¬ 
cussed so far. However, a topologically trivial gauge structure 
for density matrices has been discovered by Uhlmann [16-18] 
as we will review now. Later on (see Section III C), a more re¬ 
strictive gauge structure that can be topologically non-trivial 
will be derived from Uhlmann’s general construction. 


C. Parallel transport and geometric phases for density 
matrices 

A density matrix p is a positive semi-definite operator with 
unit trace on a Hilbert space %. Here, we consider dirnH = 
N < 00. For the application to gapped band structures that 
we have in mind here, N plays the role of the total number of 
bands of a model system, p can be represented as a pure state 
\ip) on an extended Hilbert space 7La®Rb with 7~La — 7~Lb — 
7~L such that Tr [O <8> l|'t/>)('0|] = Tr^ [Op] for any operator O 
acting in 77, a prescription referred to as state purification. 
Equivalently, the purification can be represented in terms of a 
Hilbert Schmidt operator, i.e., a N x N matrix w that satisfies 

P = Trs [|V’)(V’|] = ww\ (9) 

where the trace over the auxiliary Hilbert space 7-Lb becomes a 
matrix multiplication. This description contains a redundancy 
since under w —> wU with U G U(N ), p —> wUU^w^ = p 
is unchanged. We note a formal analogy to the basis indepen¬ 
dence of the projection ( 6 ). Furthermore, the inner product 
(w, v) = Tr[-uAi>] for matrices again (cf. Section II A) defines 
a natural means to measure the length L = J ( dr yj(w, w) 
of a path t 1 —^ w(t) of Hilbert Schmidt operators that purify a 
path t 1 —^ pit) in the sense of Eq. (9). The shortest possible 
path can again be obtained by decomposing the tangent vec¬ 
tor w into a vertical part wy and a horizontal part ib u with 
the help of the inner product, w is purely horizontal for the 
shortest path. Vertical vectors can be defined in terms of Eq. 
(9) as tangent vectors to curves s 1 —> w{s) that project to the 
same p = wvr for all s. Vertical vectors are hence of the 
form w v = jsWUis )\ s _ 0 = wu with U(s) = e“ s G U(N). 
Horizontal vectors wh are defined as orthogonal to all ver¬ 
tical vectors, i.e., Trfw^wM] = 0 for all U{N) -generators 
u = —u'. Adding the hermitian conjugate to this condition 
yields Tr[u(w* H w — w^wh)] = 0. This can be true for every 
antihermitian u only if the second antihermitian factor is zero 
by itself, i.e., the tangent vector w is horizontal if [16-18], 

w'w — w^w = 0 . 

Eq. (1), repeated here for convenience, is the density matrix 
analog of Eq. (4). 
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Let us now construct the (./(iV)-gauge field associated with 
this parallel transport prescription, in a form that is amenable 
to practical calculations. For a given p = YhjPi |j)0’|> the 
uniquely defined square root yfp = ^ |j)(i defines a 

generic purification w = y/p. In fact every purification can by 
means of a polar decomposition be written as 

w = yJ-pXJ with U G U(N ) (10) 

in some analogy to the polar decomposition of a complex 
number. Let us consider a loop 7 : t 1 —>■ p(t), t £ [0,T] 
in the space of density matrices. We denote by t 1 —>• wit) = 
yjp{t)U (f) the parallel-transport of an arbitrary initial purifi¬ 
cation w(0 ) = yjp(t))U (0) along this loop 7 and call the geo¬ 
metric phase Hijj = U {T)U{ 0)1 its Uhlmann holonomy. With 
t i = iT / M , Mg N, we can express as 

M 

Hl= lim TT U{ti)U\ti-i), ( 11 ) 

M—>00 

2=1 

where the product is ordered from right to left with increasing 
i. We now wish to obtain the analog of Eq. (7) in an explicit 
form. To this end, we first obtain the parallel-transported U ( t ) 
in Eq. (10) at infinitesimally neighbouring points in time. Eq. 
( 1 ) can be shown to be equivalent to 

w\t + e)w(t) = w\t)w(t + e) > 0 ( 12 ) 

to leading order in e; in particular, the hermitian matrix 
w' (t + e)w(t) is positive semi-definite. Defining V = w'(t + 
e)w(t) = W{t + e)y/p(t + e) yjp(fjU (t) > 0 we immedi¬ 
ately verify 

y/ p(t + e) v / p(i) = HU(t + e)U\t ) 

with the semi-positive hermitian matrix H = U(t + 
e)V ( t)U t it + e). Hence, the singular value decomposition 

\J pit + e)^/ pit) = LDR ] 

with unitary L , R yields II = LI) L' and, more importantly 
[24], 

U(t + e)(7 t (f) = LR). (13) 

With this recipe the connection Ajj = —U(i)W(t) = 
lim^o ^{l — U(t + e)W(t)) can be computed explicitly. The 
desired holonomy then reads as 

H{j = T e~ J: > Au . (14) 

We note that Au is invariant under a rescaling pit) -P 
\(t)p(t) with a strictly positive real function A(f) > 0, which 
affects the singular values D alone. Therefore, the normaliza¬ 
tion of the density matrix is not relevant for geometrical and 
topological considerations. 


III. TOPOLOGICAL ASPECTS OF MIXED STATES 

In order to discuss topological features of emergent gauge 
theories for mixed states, it is natural to view the above de¬ 
composition of tangent vectors into vertical and horizontal 
components as a connection on a principle fiber bundle (PFB) 
(see e.g.. Ref. [32]). 

It is well known that the PFB of general Hilbert-Schmidt 
operators is topologically trivial since p 1 —>■ yfp is a global 
section (gauge). At first sight, this may seem discouraging. 
On the other hand, even the special case of pure states can be 
viewed as a topologically trivial gauge structure in terms of 
Hilbert Schmidt operators. In contrast, it is well known that 
the gauge structure for pure states discussed in Section IIA 
can very well be topologically non-trivial, phenomena such 
as topological insulators being striking ramifications of this 
possibility. More concretely, the frame bundle of orthonormal 
frames ff j spanning the projection P(R) can be topologically 
non-trivial. However, it can always be viewed as a sub-bundle 
of the topologically trivial Hilbert-Schmidt bundle of arbitrary 
quadratic matrices w(R) that satisfy P(R) = w(R)w 1< (R). In 
other words, the level of description determines whether topo¬ 
logical aspects of pure states can be resolved. More generally 
speaking, every gauge theory (PFB) can be viewed as a sub¬ 
bundle of a topologically trivial gauge theory (PFB) with a 
bigger gauge degree of freedom but over the same base man¬ 
ifold. As for mixed state density matrices key questions are 
thus whether topologically non-trivial features exist and, if so, 
to identify a gauge structure which is able to reveal their topo¬ 
logical content. Here we consider several physically moti¬ 
vated constraints under which a topologically non-trivial PFB 
of EOPS will be defined. 


A. Triviality and structure of the Hilbert-Schmidt bundle 

Let us first review the case of arbitrary invertible N x N 
density matrices p as considered by Uhlmann [16], i.e., strictly 
positive matrices without any further constraints such as nor¬ 
malization of the trace. The space Mm{N) of all Hilbert 
Schmidt operators w purifying such density matrices via Eq. 
(9) is then simply given by the group of all invertible matri¬ 
ces GL(C,N). This is because p = ww^ is strictly posi¬ 
tive and has the non-vanishing determinant |det(w )| 2 if and 
only if w G GL(C,N). Mn(N) = GL(C,N) is the to¬ 
tal space of the Hilbert-Schmidt bundle which projects via 
H : w 1 —^ ww' onto the base manifold of strictly positive 
matrices Dn(N) = GL(C,N)/U(N). The quotient form of 
Dn{N) is rooted in the fact that any regular matrix M can 
be uniquely written as M = \JMM t U. U G U(N). For the 
same reason, w with ww' = p, can be uniquely represented 
as w = yJpU, U G U(N). Hence the fiber over p is given by 
G p = { y/pU : U G U(N)J, which is manifestly isomorphic 

to U(N). GL(C, N) Dn(N) thus defines a PFB. Due to 
the existence of the global gauge or section p —>■ yfp, the PFB 
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is topologically trivial, i.e., 

GL(C,N) = GL{C,N)/U{N) x U(N). (15) 

Now we consider singular density matrices, i.e., semi-positive 
N x N matrices with rank n < N. A naive analog of the 
regular case with fibers G p = {y/pU : U £ U(N)} where p 
now has rank n does not give a PFB. The reason is that the 
right-action w = ^fpU of U(N) is no longer free, i.e., many 
U £ U(N) give the same w = yfpU due to the N — 71- 
dimensional null-space of ^fp. Simply restricting the U to 
U (n) acting on the support of ^fp does not resolve this issue 
since this action would not be transitive, i.e., unable to reach 
all w with ww' = p. In Ref. [33], it has been shown that 

Mjv(n) -5- n < N with the same projection II : 

w 1 —f wv/' still defines a PFB. However, by the same argument 
as before, namely by the existence of the global section p 1-4 
y fp, also these bundles are topologically trivial. 

B. Ensemble of pure states gauge structure 

Instead of considering general Hilbert Schmidt operators, 
we define a EOPS bundle in terms of n non-orthonormalized 
pure states \ip a ) that project onto a given N x N density ma¬ 
trix p of rank n as (Eq. (2) is repeated for convenience) 

n 

p = £ir)<ri 

a=l 

as a gauge structure over the space of density matrices p. In¬ 
stead of viewing p as a single pure state in a larger Hilbert 
space, EOPS in this scheme are plausible ensembles of pure 
states in the system Hilbert space. A generic choice | i/j a ) = 
y/p^\ip a ) along these lines is obtained from the spectral rep¬ 
resentation p = '^2 a Pa\'4 , a){4 > a\- Under arbitrary unitary ro¬ 
tations 

I ‘Pa) = |V>/3 )Up a , (16) 

the density matrix p obtained from the projection ( 2 ) is un¬ 
changed. However, the pure states \p a ) are no longer mutu¬ 
ally orthogonal. The matrix representation 

«j=(|^i),...,|^ n ),0,...,0) (17) 

with N — n zero columns defines a natural embedding of 
the space Pm (n) of non-orthonormal n-frames into the space 
Mjv(n) of Hilbert Schmidt operators of rank n [33], U as 
occurring in Eq. (16) defines the local U(n ) gauge degree 

of freedom of the PFB Pjv(n) Dm(ti). For regular den¬ 
sity matrices (n = N ), it is easy to see that the space of all 
EOPS is identical to the space of all Hilbert Schmidt oper¬ 
ators satisfying p = wwK i.e., Pm{N) = Mn(N). Thus, 
for invertible density matrices with no further spectral con¬ 
straints, the EOPS scheme is equivalent to the purification by 
pure state representation in a larger Hilbert space. From this 


observation, we immediately conclude that the EOPS bundle 
is topologically trivial for unconstrained regular density ma¬ 
trices. In the following, we will consider several constraints 
under which the EOPS bundle can become topologically non¬ 
trivial and explicitly construct the corresponding topological 
invariants. 


C. Gauge structure of spectrally constrained density matrices 

Non-degenerate density matrices - The situation becomes 
more interesting if we impose certain conditions on the spec¬ 
trum of the density matrix. Let us first consider a regular non¬ 
degenerate density matrix p = )('*/’«I- where we 

can now without loss of generality assume pi > p 2 > ■.. > 
Pm- For the hermitian operator p with non-denerate eigenval¬ 
ues, the projections P a = jpa) (Vt*! are mutually orthogonal 
and it is natural to order the columns y/p^lPa) of a Hilbert 
Schmidt representation (17) with descending size of the spec¬ 
tral weights p a . This ordering will only be maintained under 
a subgroup of gauge transformations (16) that are a direct sum 
of N (7(1)-transformations e 1 ^" acting on the rays of eigen¬ 
states associated with the eigenvalues p a . Eq. (16) then sim¬ 
plifies to 

\(p a ) = a = l,...,N. (18) 

The EOPS bundle thus constrained consists of N (7(1)- 
bundles Pjy(l) ^ ipv(l) which are subject to the constraint 

N 

y^lV , o;)(^a| = 1- (19) 

a—1 

While the individual (7(l)-bundles can be topologically non¬ 
trivial, Eq. (19) enforces a zero sum rule for their topolog¬ 
ical invariants. The reason for this is analogous to the pure 
state case and can be intuitively understood as follows. The 
individual subspaces may exhibit a topologically non-trivial 
winding as a function of the control parameters (e.g. lattice 
momentum) in the total Hilbert space. This total space serves 
as fixed reference and thus by definition does not exhibit any 
’’net-winding”. Hence, the parameter dependences of the in¬ 
dividual subspaces have to compensate each other in order to 
span the resolution of the identity (19) of the total embedding 
space at every point in parameter space. 

The spectral constraint forbidding any level degeneracy can 
be somewhat relaxed by assuming that m subsets of cardinal¬ 
ity ni, ..., n m with N = X]j=i n j °f levels are allowed to 
be degenerate, but still have distinct eigenvalues from levels 
outside of their subset. Such constraints simply result in a 
restricted gauge group of the form U(n\) x U(ri2) x ... x 
U(n m ). Again, the individual U(rij) -structures may be topo¬ 
logically non-trivial but obey a zero sum rule due to Eq. (19). 

Purity gaps - The notion of a purity gap [21, 22] can be 
rationalized in this framework. To this end, we consider a 
spectral constraint along the lines of the above discussion 
where the n largest eigenvalues of p belong to one subset 
but are by assumption not degenerate with the (n + l)-th 
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largest eigenvalue of p. Under these circumstances we can 
uniquely define a pure state projection P = X^a=il'0a)('0a| 
which is closest to the mixed state defined by p. In Refs. 
[ 21 , 22 ], the topological invariant of p has been defined as 
the invariant of the corresponding pure state projection P. 
Below, we directly define topological invariants for mixed 
states in terms of the Uhlmann connection (12) restricted to 
the subspace associated with the n largest eigenvalues of p. 
Their value is equal to the one obtained by performing an 
adiabatic deformation into a pure state along the lines of Refs. 
[ 21 , 22 ] but, quite remarkably, no reference to pure states is 
required here. A purity gap closing means that the n-th largest 
eigenvalue and the [n + l)-th largest eigenvalue of p become 
degenerate, a singular-point at which the spectral constraint 
is violated and the topological invariant associated with the 
subspace of the n largest eigenvalues is not well defined. 
As a consequence, this subspace can exchange topological 
charge with the eigenspace associated to the (n+ l)-th largest 
eigenvalue and the topological invariant for the subspace of 
the n largest eigenvalues may have changed when the purity 
gap reopens - a topological phase transition by purity gap 
closing. 


of interest by virtue of V. We denote by llj/"' the repre¬ 
sentation matrix of the geometric phase in this subspace, i.e., 

^ H\ipj). With these definitions, the cur¬ 
vature of the U ( n ) EOPS bundle is given by 


P 




lim 

a—>o 


i - &ir v 


( 20 ) 


where <5 —0 is shorthand for (<5^, 5 V ) (0, 0). In particular, 

the gauge invariant Abelian curvature reads as 


■BA T f v logdet {H% tv ) 

P^ = TrJ^ = - lim 


<5—>o 


Sy.6 u 


( 21 ) 


where logdet(.H r ^j“') is i times the argument of the Abelian 
geometric phase factor. In terms of this curvature, stan¬ 
dard topological invariants for mixed states with spectral con¬ 
straints can be defined. For example, if the constrained rank n 
density matrix p(k) is parameterized by a lattice momentum 
k in a two-dimensional Brillouin zone (BZ), the first Chern 
number is given by 


Singular density matrices - A physical system under non¬ 
equilibrium conditions may not be ergodic and can thus be 
described by a singular density matrix, i.e., a density matrix 
of rank n < N. This constraint can be viewed as special case 
of the above discussion of spectral constraints: A subset of 
n eigenvalues p a is non-zero and hence not degenerate with 
the complement of N — n non-zero eigenvalues. The U(n) 
EOPS bundle associated with the non-zero subspace can be 
topologically non-trivial. 


D. Topological invariants 

The gauge structure of spectrally constrained density matri¬ 
ces defined in Section III C allows for the definition of topo¬ 
logical invariants directly in terms of the Uhlmann connection 
(see Section II C), i.e., without reference to pure states. To 
this end, we consider a subset of n < N eigenvalues (eigen¬ 
vectors) labeled by j = 1 ,,n that are by assumption non¬ 
degenerate with (orthogonal to) the N—n remaining eigenval¬ 
ues (eigenvectors) of p. We denote by V = ^" =1 I V’j) {^j I the 
projection onto the rank n density matrix p = JA1 tpj)('tpj\, 
i.e., p = VffP. Eq. (17) shows how the corresponding 
U(n) EOPS bundle can be embedded into the general PFB 
of Hilbert Schmidt operators considered by Uhlmann. This 
allows us to construct a curvature on the constrained U(n ) 
EOPS bundle which is directly inherited from Uhlmann’s gen¬ 
eral construction. By general definition, a curvature P /n/ is the 
geometric phase per area associated with the parallel trans¬ 
port around and infinitesimal parallelogram y lu/ spanned by 
the vectors in momentum space (or more gener¬ 

ally in parameter space). To obtain the curvature of the con¬ 
strained EOPS bundle, we have to project Uhlmann’s geo¬ 
metric phase (see Eq. (14)) associated with the in¬ 

finitesimal loop 7 ;|I/ for p onto the n-dimensional subspace 


C = 




iP A 

l J xy 

27T 


( 22 ) 


The generalization to other invariants and the implementation 
of physical symmetries is straight forward and analogous to 
the pure state case. As a prominent example, we would like 
to mention the case of time reversal symmetry as a protect¬ 
ing symmetry (see Section IIB for a discussion of the pure 
state analog). Instead of the projection onto the occupied 
bands, the density matrix p{k) itself then obeys the symme¬ 
try Tp{k)T ~ 1 = p{—k) with the anti-unitary time reversal 
operator T. As in the pure state case, this symmetry carries 
over to arbitrary Uhlmann-holonomies and allows for the def¬ 
inition of topological invariants in analogy to Ref. [34], 


IV. UHLMANN PHASE WINDING NUMBERS 

In two recent back-to-back publications [24, 25], a com¬ 
plementary approach towards the definition of topological in¬ 
variants for mixed states has been reported, in particular for 
Chern insulators in thermal states. Their approach is based on 
the so called Uhlmann phase (b[j associated with a loop 7 in 
momentum space, defined as 

= Tr [w(0) t rt;(T)] = Tr [p(0)U(T)U{ 0) f ] 

= Tr \p(0)Hy] , (23) 

where 7 is traversed between t = 0 and t = T, w(t) = 
yjp(t)U(t) is a parallel transport in the sense of Uhlmann’s 
connection (1) and the Uhlmann holonomy H' l - has been de¬ 
fined in Eq. (14). By construction, (fy bears some anal¬ 
ogy with the Abelian Berry phase of a pure state defined 
as e 1 ^ = ('i/'(0)|-i/>(T)) (cf. Section II A). However, there 
is one crucial difference in the mathematical stmcture of © 7 
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and (jyjj-. While t 1 ^ is a [/(l)-holonomy as is clear from the 
explicit representation in terms of an integral over a U( 1)- 
gauge field (see Eq. (5)), does in general not have 
this property - in contrast to the full non-Abelian Uhlmann 
holonomy H^. This observation is not just a minor techni¬ 
cal point but can have drastic consequences. The key fea¬ 
ture of a (7(1)-holonomy is its additive group structure, i.e., 
cA 712 = </> 71 +(P 2 (mod27r), where 712 denotes the concatena¬ 
tion of the loops 71 and 72 . It is this very property that allows 
us to go from an infinitesimal geometric phase as measured 
locally by a curvature (see Eq. (20)) to a quantized global 
topological invariant represented as an integral over the entire 
parameter space (see, e.g., the definition of the Chern num¬ 
ber in Eq. (22)). Again making use of this additive group 
structure of the ordinary Berry phase, the Chern number C 
can be represented as the change of the Berry phase (!> k:r = 
f gl dk y A y (k x , ky) associated with a 7’ y -circle in the BZ at 
fixed k x . This reformulation is achieved by dividing the BZ 
torus into fc y -rings of infinitesimal width [k x , k x + dk x ] 
around k x and then using the Stokes theorem to transform 
the area integral of the curvature over these rings into a 
line integral of the connection A along their boundaries, i.e., 

,[/j. tl kj~ r y J S | (iky fTlyy ( k X T (lk X , ky ) A.y (/>7 ) (*'7/ ) ] 

(t»T~) ( ^ x ' Explicitly, this procedure leads to the following 
well known representation of the Chern number 



(24) 


where the circle S' 1 is around the k x loop of the BZ. Equiv¬ 
alently, the role of k x and k y could be exchanged here, of 
course. This is again because of the holonomy group struc¬ 
ture of the Berry phase which allows us to divide the BZ into 
stripes in an arbitrary direction. 

The basic idea of Refs. [24, 25] is to take the right hand 
side of Eq. (24) but to replace the ordinary Berry phase by the 
Uhlmann phase as defined in Eq. (23), i.e., 


and k y in all calculations, i.e. 



FIG. 1. (color online) Momentum dependence of the Uhlmann phase 
for a thermal state at /3 = 1.3 for the model defined in Eq. (27), 
numerically obtained by discretizing the paths in momentum space 
into 500 equidistant points. Left panel: cj> L f as a function of k x . 
k i 

Right panel: (j)J as a function of k y . 


Here, we demonstrate that indeed Cu A Cjj can occur for 
a Chern insulator the Hamiltonian of which is not invariant 
under the exchange of k x and k y . To this end, we consider a 
two-banded fermionic Chern insulator on a 2D square lattice 
with unit lattice constant defined by the Bloch Hamiltonian 

3 

H(k ) = d(k) • a = d? (k)aj, 

i 

d 1 (k) = sin (k x ), d 2 (k) = 3 sin (k y ), (27) 

d 3 {k) = 1 — cos(fca;) — cos (k y ), 

where cr, are Pauli matrices. Note the anisotropy factor of 
3 in d 2 (k) which is crucial here. For the two-banded model 
(27), the Uhlmann connection and the Uhlmann phase can be 
readily calculated along the lines of Ref. [35]. Explicitly, the 
Uhlmann connection An in this special case reads as 


Cu = 



d$j_ 

dk x 


d k T . 


(25) 


We note that Ref. [24] actually discusses several constructions 
in terms of the spectrum of the “holonomy matrix” p(0)H^. 
Since the absence of a holonomy group structure which is at 
the heart of our present discussion also pertains to p(0)H^ r , 
this distinction is not of central importance here as will be 
addressed more explicitly in our example below. As both k x in 
the BZ and the phase (j^f are defined on a circle, the mapping 
k x i—^ <\njj is characterized by an integer quantized winding 
number which is exactly measured by Eq. (25). 

However, since the Uhlmann phase does not have an addi¬ 
tive group structure, Cu cannot be represented as the integral 
of a curvature over the 2D BZ. It is hence not immediately 
clear to what extent Cu, technically being a ID winding num¬ 
ber, can be seen as a unique property of the 2D system under 
investigation. In particular, it is not clear that Cu as defined in 
Eq. (25) is equal to Cu obtained from Cu by exchanging k x 


Au,y, 


(^t/#)). \fpW> 


(28) 


At T = 0, Eq. (28) concurs with the Kato connection [13] 
(see Eq. ( 8 )) for the pure state case and the ordinary Chern 
number (24) reads as 

c = ^ j *k ( d{k) • x (d k j(k))]) 

with d = t 4. Explicit calculation yields C = — 1. Also, both 

~ 1^1 

Cu and Cu trivially concur with C since the Uhlmann phase 
reduces to the Berry phase for pure states such that Eq. (24), 
Eq. (25), and Eq. (26) become equivalent. At finite T, the 
system is in a thermal state defined by p(k) = -7 e -7#(fc) 
with Z = Tr[e“^ ff ( fc )]. Using the simplified form (28) of 
the Uhlmann connection, the direct calculation of Uhlmann 
holonomies (14) for arbitrary loops in momentum space is 
straightforward which in turn allows the direct calculation 
of Cu and Cu from Eq. (25) and Eq. (26), respectively. 


















8 


In the infinite temperature limit, this calculation becomes 
trivial as p(k) = 1 independent of k and hence Au = 0 in 
Eq. (28), resulting in Cjj = Cjj = 0. Hence, the Uhlmann 
phase winding numbers have to jump from —1 to 0 at some 
temperature. Numerically performing this calculation for 
finite temperature, we find that both Cjj and Cjj jump to zero 
at finite critical temperatures /3 C = jr and /3 C = i, respec¬ 
tively. Remarkably, for the present model, 8c ^ B c , i.e., 
Cjj and Cjj do not concur at all temperatures. Numerically, 
we find j3c = 0.874 and /3 c = 1.32. This discrepancy is 
illustrated In Fig. 1, where we visualize the winding numbers 
Cjj (left panel) and Cjj (right panel) for 8 = -^ = 1-3, i.e., 

8c < 8 < L Clearly, Cjj = 1 since <f>u monotonously 
increases, interrupted by a jump from 7r to — 7r at k = 0, 

by 27 t as k x completes the loop from k x = —7r to k x = n. 

k 

In contrast, cbj does not reach all values in [—7r, 7r] which 
implies Cjj = 0. On a more detailed note, we would like 
to point out that, for the parameters chosen here, Cjj is in 
the well defined regime in the sense of Ref. [24], This 
is because the “holonomy matrix” p(k x , k y o)H^j k:c ' ky0 of 
a loop with footpoint k y o at fixed k x is gapped for all 
kx,k y o, i.e., the difference of the absolute values of its 
eigenvalues is finite and bounded from below by A = 0.268. 
Comparing the construction in Refs. [24, 25] to our present 
analysis in Section III, we would like to emphasize that the 
Chern number defined through Eq. (20) and Eq. (22) does 
not exhibit similar finite temperature transitions as Cjj and 
Cjj but can, for the model (28), be uniquely defined as — 1 
for the larger eigenvalue of the density matrix for all finite 
temperatures. At infinite temperature, in contrast, the density 
matrix does not obey any spectral constraints in the sense 
discussed in Section IIIC, rendering all topological invariants 
trivial. 


V. CONCLUDING REMARKS 

We have discussed how several assumptions regarding 
the spectrum of a family of density matrices can lead to a 
topologically non-trivial gauge structure. In this framework 
topological invariants that are protected by these spectral 
assumptions have been defined for mixed states. Protected 
here means that topologically inequivalent mappings from a 
parameter space into the density matrices can be continuously 
deformed into each other only if the underlying spectral 
assumptions are violated. Identification of the parameter 
space with the Brillouin zone of a lattice translation invariant 
system provides one way of generalizing topological band 
structure invariants to the realm of mixed states. A non-trivial 
example where going beyond pure states is crucial to obtain 
a topologically non-trivial Chern insulator as a steady state 
in the framework of a non-equilibrium open quantum system 
dynamics has been reported in Ref. [36]. Additional physical 
symmetries refining this system of topological invariants 
can be considered in analogy to the pure state case. The 
topological invariants defined here, being gauge invariant 
properties of the density matrix, are in principle experimen¬ 
tally accessible via state tomography. However, their relation 
to natural observables such as response functions is not yet 
conclusively understood. First progress along these lines 
has been reported in Ref. [37]. As for two-banded Chern 
insulators such as the toy model in Eq. (27), any statistical 
mixture of bands with opposite Chern number will certainly 
cause deviations from the quantized Hall conductance. 
In other situations where the physical ramification of the 
topological invariant is related to a half-integer quantized po¬ 
larization, a statistical mixture may still exhibit a quantization. 
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